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Abstract. We present self-consistent general relativistic simulations of stellar core collapse, 
bounce, and postbounce evolution for 13, 15, and 20 solar mass progenitors in spherical symmetry. 
Our simulations implement three-flavor Boltzmann neutrino transport and standard nuclear physics. 
The results are compared to our corresponding simulations with Newtonian hydrodynamics and 
0(v/c) Boltzmann transport. 

INTRODUCTION 

A supernova explosion is a dramatic event that includes such a rich diversity of physics 
that, with current computer hardware, a self-consistent model based on numerical simu- 
lations can not possibly include all of them at once. After stellar core collapse, a compact 
object at the center of the event is formed, requiring a description in general relativity. 
Neutrinos radiating from this central object are strongly coupled to the matter at high 
densities before streaming out at lower densities. Multi-frequency radiation hydrody- 
namics must be used to quantify the energy that these neutrinos deposit in the material 
behind the shock. Moreover, evidence suggests that this heating drives convection be- 
hind the shock, and significant rotation and strong magnetic fields might also be present. 
Observations of neutron star kicks, mixing of species, inhomogeneous ejecta, and po- 
larization of spectra support the presence of asymmetries in supernova explosions (e.g. 
Tueller et al. [1], Strom et al. [2], Galama et al. [3], Leonard et al. [4], and references 
therein). Motivated by such observations, various multi-dimensional explosion mecha- 
nisms have been explored (Herant et al.[5]. Miller et al.[6], Herant et al.[7]. Burrows et 
al.[8], Janka and Miiller [9], Mezzacappa et al.[10]. Fryer [11], Fryer and Heger [12]) 
and jet-based explosion scenarios have recently received new momentum (Hoflich et al. 
[13], Khokhlov et al. [14], MacFadyen and Woosley [15], Wheeler et al. [16]). Many 
exciting aspects of jets are discussed in these proceedings. However, for the time being, 
one has to single out a subset of the known physics and to investigate the role each part 
plays in a restricted simulation. It is natural to start with ingredients that have long been 
believed to be essential for the explosion and to add modifiers in a systematic way until 
the observables can be reproduced (see e.g. Mezzacappa et al. [17]). Spherically sym- 
metric supernova modeling has a long tradition and is nearing a definitive point in the 



sense that high resolution hydrodynamics, general relativity, complete Boltzmann neu- 
trino transport, and reasonable nuclear and weak interaction physics are being combined 
to dispel remaining uncertainties. 

Our simulations are based on the three-flavor Boltzmann solver of Mezzacappa and 
Bruenn [18, 19, 20] that was consistently coupled to hydrodynamics on an adaptive 
grid and extended to general relativistic flows (Liebendorfer [21]). We performed rela- 
tivistic simulations for progenitors with different masses and compare the results to the 
counterpart simulations in Newtonian gravity with 0(v/c) Boltzmann transport (Messer 
[22]). While the quantitative results are sensitive to the inclusion of all neutrino flavors 
and general relativity (Bruenn et al. [23]), we find the same qualitative outcome as did 
Rampp and Janka [24] in their independent simulations with one-flavor 0(v/c) Boltz- 
mann transport: Given current input physics, there are no explosions in spherical sym- 
metry without invoking multidimensional effects (Mezzacappa et al. [17], Liebendorfer 
et al. [25]). This is consistent with the results of Wilson and Mayle [26], who found 
explosions only by including neutron finger convection. 



CORE COLLAPSE AND NEUTRINO BURST 

The simulations discussed herein are initiated from progenitors with main sequence 
masses of 13 Mq, 20 M0 (Nomoto and Hashimoto [27]), and 15 M0 (Woosley and 
Weaver [28]). We use in our models the equation of state of Lattimer and Swesty [29] 
and "standard" weak interactions (e.g. Bruenn [30]). The first phase in the simulations 
leads through core collapse to bounce and shock formation. When the shock passes 
the neutrinospheres approximatively 4 ms after bounce, an energetic neutrino burst is 
released from the hot shocked material, rendering it "neutrino-visible" to the outside 
world. The corresponding deleptonization behind the shock is dramatic. The energy 
carried off with the neutrino burst enervates the shock in both the NR and the GR 
cases. A pure accretion shock continues to propagate outwards as infalling material is 
dissociated and layered on the hot mantle. This stage is the definitive end of a "prompt", 
i.e. purely hydrodynamic, explosion. 

The neutrino luminosities from the general relativistic simulations are shown in Fig. 
(1). The electron neutrino luminosity slowly rises during collapse and decreases as the 
core reaches maximal density. It remains suppressed for the ~ 4 ms the shock needs 
to propagate to the electron neutrinosphere. The most prominent feature is the electron 
neutrino burst, reaching 3.5 x 10^^ erg/s at shock breakout, and declining afterwards. 

The initial similarities in density and temperature in the inner parts of the three cores 
are responsible for the similar evolution of each core with respect to collapse, bounce, 
and the signature of the electron neutrino burst (Messer [22]). This is shown in Figs. 
(1) and (2). Any initial difference in the electron fraction, Yg, in the inner part of each 
core cannot survive in regions that reach sufficient densities for neutrino trapping and 
equilibration to occur. If a mass element reaches neutrino trapping density, the total 
lepton fraction, 7;, becomes the variable determining the state of the mass element. At 
equilibrium, the electron capture and neutrino absorption rates are related to one another 
through detailed balance, and the final Yg becomes a function of the local temperature 




FIGURE 1. Luminosities for all flavors in the GR simulations. The last graph shows an energy rate 
calculated by the product of the gravitational potential at the neutrinosphere and the mass flux crossing it. 
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FIGURE 2. Similar core profiles in the GR cases at 10 ms after bounce. 



and density through the electron chemical potential. Since the inner temperature and 
density profiles differ very little between the models, each settles to essentially the same 
equilibrium 7^ deep in the core. This determines the size of the homologous inner core 
and, consequently, the position of shock formation. Since this position sets the amount 
of mass that has to be dissociated when the shock ploughs to the neutrinospheres, it 
strongly influences the shock energy. Thus, a strikingly similar shock initiation is found 
for the three progenitors. However, the initially dissociated mass differs largely between 
the GR and NR cases because of the smaller enclosed mass at shock formation due to the 
GR effects in the gravitational potential (GR: 0.53 Mq, NR: 0.65 Mq). The luminosities 
in the GR and NR cases are very similar for the electron neutrino burst. Afterwards, the 
GR luminosities are generally 10% — 20% larger than the NR luminosities. 



HEATING AND ACCRETION 

In the standard picture, the ensuing evolution is driven by electron neutrino heating. 
Electron flavor neutrinos diffuse out of the cold unshocked core and are created in the 
accreting and compactifying matter around the neutrinospheres in a hot shocked mantle. 
If the heating is not sufficient to revive the shock within the first ~ 200 ms after bounce 
(as is the case in our simulations), the material in the heating region is drained from 
below (Janka [31]), and the conditions for heating deteriorate because of the shortened 
time the infalling matter spends in the heating region. We observe this process earlier in 
the general relativistic simulations than in the nonrelativistic simulations. 

The shock trajectories for the different progenitors are shown in Fig. (3). The main 
difference between the GR and NR simulations stems from the difference in the size of 
the proto-neutron star, which is caused by the nonlinear GR effects at very high densities 
in the center of the star. At radii of order 100 km and larger, the gravitational potential 
becomes comparable in the GR and NR cases. However, large differences arise if the 
steep gravitational well is probed at different positions in the GR and NR simulations, 
as happens with accretion down to the neutrinospheres. The deeper neutrinospheres 



300 



250 



200 



Shock Trajectories 



E 



150 



T3 

to 



100 



50 







GR, 13Ms, Nomoto and Hashimoto 88 

GR, 15Ms, Woosley and Weaver 95 

GR, 20Ms, Nomoto and Hashimoto 88 

NR, 13Ms, Nomoto and Hashimoto 88 

NR, 15Ms, Woosley and Weaver 95 

NR, 20Ms, Nomoto and Hashimoto 88 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Time After Bounce [s] 



FIGURE 3. Shock trajectories for the GR and NR cases. 



in the more compact GR case encounter material that has traversed a larger potential 
difference, producing more energetic accretion luminosities. After the shock has stalled, 
the surrounding layers (cooling/heating region, shock radius) adjust to a smaller radius 
and settle to a stationary equilibrium in the spirit of Burrows and Goshy [32]. In the GR 
simulation, this occurs at a smaller radius, deeper in the gravitational well, with higher 
infall velocities, leading to higher accretion luminosities with higher rms energies and 
higher heating rates (see also Bruenn et al. [23]). 

This accretion-determined picture obtains additional support from the following ob- 
servation: In the last graph in Fig. (1) we provide the mass flux at the electron neutri- 
nosphere multiplied by the local gravitational potential. We can compare this energy rate 
to the luminosities of the electron flavor neutrinos that are shown in the second graph 
in Fig. (1). The striking similarity suggests that the luminosity is indeed mainly deter- 
mined by the gravitational potential at the boundary of the proto-neutron star and the 
mass accretion rate. Moreover, we can directly relate the different accretion rates to dif- 
ferences in the density profiles of the initial progenitors. The energy accumulation rate 
at the neutrinospheres for the 15 Mq model exceeds that of the 20 Mq model between 
60 ms and 160 ms after bounce. We have marked in Fig. (2) the region that crosses the 
neutrinospheres during this time. The density of the 15 Mq progenitor exceeds the den- 
sity of the 20 Mq progenitor in exactly this region. Thus, the time dependent electron 
flavor luminosities reflect the position of the neutrinospheres and are directly modulated 
by the variations in the spatial density profiles of the progenitors. 
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